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Through the triple-alpha process occurring in red giant stars the bulk of the carbon 
existing in our universe is produced. We calculated the change of the triple-alpha 
reaction rate for slight variations of the nucleon-nucleon force using a microscopic 
12-body model. Stellar model calculations for a low-mass, intermediate-mass and 
massive star using the different triple-alpha reaction rates obtained with different 
strengths of the N-N interaction have been performed. Even with a change of 0.4% 
in the strength of N-N force, carbon-based life appears to be impossible, since all 
the stars then would produce either almost solely carbon or oxygen, but could not 
produce both elements. 



1 Introduction 

The Anthropic Principle in its weak form can be formulated in the following 
way: The observed values of all physical and cosmological quantities are not 
equally probable but they take on values restricted by the requirement that there 
exist .sites where carbon-based life can evolve and by the requirement that the 
universe be old enough for it to have already done sou. The triple-alpha process 
occurring in the helium burning of red giants is of special significance with 
respect to the anthropic principle. Through this process the bulk of carbon, 
e.g., the carbon isotope ^^C, existing in our universe is produced. Carbon is 
synthesized further by alpha capture to oxygen, e.g., to the oxygen isotope 
^^O, leading to abundance ratio of ^^C : ^^O « 1 : 2. Without the production 
of an appreciable amount of carbon, obviously no carbon-based life could have 
developed in the universe. On the other hand, the production of oxygen is also 
absolutely necessary, because no spontaneous development of carbon-based life 
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is possible without the existence of water (H2O). 

The reason for the relevance of the triple-alpha process with respect to 
the anthropic principle lies also in the fact that one has to deal with physical 
quantities that are in the realm of experimentally verifiable and theoretically 
calculable physics. This is for instance hardly the case for the much less well- 
known and complicated science necessary for the description of the Big Bang 
as well as for the creation and evolution of life on earth. 

The formation of ^^C in hydrogen burning is blocked by the absence^ 
stable elements at mass numbers A = 5 and A = 8. Opik and Salpeter ElB 
pointed out that the lifetime of ^Be is long enough, so that the a + a ^ *Be 
reaction can produce macroscopic amounts of equilibrium ^Be in stars. Then, 
the unstable ®Be could capture an additional a-particle to produce stable ^^C. 
However, this so-called triple-alpha reaction has very low rate since the density 
of ^Be in the stellar plasma is very low, because of its short lifetime of 10~^^ s. 

Hoyleu argued that in order to explain the measured abundance of carbon 
in the Universe, the triple-alpha reaction cannot produce enough carbon in a 
non-resonant way, therefore it must proceed through a hypothetical resonance 
of ^^C, thus strongly enhancing the cross section. Hoyle suggested that this 
resonance is a = 0+ state at about e = 0.4 McV (throughout this paper 
s denotes resonance energy in the center-of-mass frame relative to the three- 
alpha threshold, while F denotes the full width). Subsequent experiments 
indeed found a 0^ resonance in ^^C in the predicted £nergy region by studies of 
the reaction ^^N(d,Q!)^^C and the /ir-decay of ^^BqeI. It is the second O'^ state 
(Oj) in i^C. Its modern parametersO, e = 0.3796 MeV and T = 8.5 x 10"^ MeV, 
agree well with the old theoretical prediction. 

In this work we investigate the amount of carbon and oxygen production 
in helium burning of red giants by slightly varying the nucleon-nucleon inter- 
action. In Ref. Q only hypothetical ad hoc shifts of the resonance energy of the 
Oj-state were investigated, whereas in this work we start by variations of the 
underlying N-N interaction. 

After the Introduction we discuss in Sect. 2 our microscopic three-cluster 
model, the effective nucleon-nucleon (N-N) interactions, and the complex scal- 
ing method, which is used to describe the Oj state of ^^C. In Sect. 3 the 
change of the triple-alpha reaction rates by slightly varying the underlying 
N-N interaction is calculated. In Sect. 4 we present the results for stellar 
model calculations for a low-mass, intermediate-mass and massive star using 
the triple-alpha reaction rates obtained with different strengths of the effective 
N-N interaction. In Sect. 5 the obtained results are summarized. 
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2 Nuclear physics 

The astrophysical models that determine the amount of carbon and oxygen 
produced in red giant stars need some nuclear properties of ^^C as input pa- 
rameters. Namely, the position and width of the 0^ resonance, which almost 
solely determines the triple-alpha reaction rate, and the radiative decay width 
for the O2 2^ transition in ^^C. Here we calculate these quantities in a 
microscopic 12-body model. 

In order to make such a calculation feasible, we use the microscopic cluster 
model. This approach assumes that the wave function of certain nuclei, like 
^^C, contain, with large weight, components which describe the given nucleus 
as a compound of 2-3 clusters. By assuming rather simple structures for the 
cluster internal states, the relative motions between the clusters, which are 
the most important degrees of freedom, can be treated rigorously. The strong 
binding of the free alpha-particle CHel makes it natural that the low-lying 
states of ^■^C are largely 3a-structures Therefore, our cluster-model wave 
function for ^^C looks like 

*"^-E-^{'J'"*"*"^M.!i(Pi.P2)}- (1) 

Here A is the intercluster antisymmetrizer, the cluster internal states are 
translationally invariant Os harmonic-oscillator shell-model states with zero 
total spin, the p vectors are the intercluster Jacobi coordinates, li and I2 are 
the angular momenta of the two relative motions, L is the total orbital angular 
momentum and [. . .] denotes angular momentum coupling. The total spin and 
parity of ^^C are J = L and tt = (— respectively. 

We want to calculate the resonance energy of the state in ^^C, relative to 
the 3a-threshold, and the 0^ 2^ radiative (E2) width, while slightly varying 
the strength of the effective nucleon-nucleon (N-N) interaction. This way, the 
sensitivity of the carbon and oxygen production on the effective N-N interaction 
can be studied. The 0^ state is situated above the 3a-threshold, therefore for a 
rigorous description one has to use an approach which can describe three-body 
resonances correctly. We choose the complex scaling methodEl that has already 
been used in a variety of other nuclear physics problems, see e.g. cB. In this 
method the eigenvalue problem of a transformed Hamiltonian 

He = u{e)HU-\e) (2) 

is solved, instead of the original H , where the transformation U acts on a 
hmction /(r) as 

U{e)f{v) = e3*«/2/(re*'). (3) 
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Resonances are eigenvalues of H with i^ics = e — iV /2 (e, F > 0) complex 
energy, where e is the position of the resonance while T is the full width. A 
straightforward description of such states are difficult because their wave func- 
tions are exponentially divergent in the asymptotic region. The effect of the 
Eq. complex scaling transformation is that the positive-energy continuum 
of H gets rotated (by the angle 20) down into the complex energy plane, while 
the wave function of any resonance becomes square-integrable if 26* > | arg E^-cs \ ■ 
Using this methodtl, we were able to localize the energy of the Qt state of ^^c. 

The other important quantity that needs to be calculated is the radiative 
width of the 0^ state, coming from the electric dipole (E2) decay into the 
2'^ state of ^^C. This calculation involves the evaluation of the E2 operator 
between the initial J three-body scattering state and the final 2^ bound state 
Eil. The proper three-body scattering-state treatment of the initial state is 
not feasible for the time being, therefore we use a bound-state approximation 
to it. This is an excellent approximation for the calculation of F-y because the 
total width of the state is very small (8.5 eVQ). The value of F^ is rather 
sensitive to the energy difference between the 0^ and 2^ states, so we have to 
make sure that the experimental energy difference is correctly reproduced. 

In order to see the dependence of the results on the chosen effective N-N 
interaction, we performed the calculations using three different forces. The 
Minnesota (MN)ll3, and the rather different Volkov 1 and 2 (VI, ¥2)113 forces 
achieve similar quality in describing light nuclear systems. We slightly adjusted 
a parameter (the exchange mixture) of each forces in order to get the 
resonance energy right at the experimental value. This leads to m = 0.941 for 
MN and m — 0.568 and m — 0.594 for VI and V2, respectively. The results 
coming from these forces are our baseline predictions. Then, we multiplied 
the strengths of the forces by a factor p, which took up the values p = 0.996, 
0.998, 0.999, 1.001, 1.002, and 1.004, respectively, and calculated the resonance 
energies and gamma widths again. This way we can monitor the sensitivity 
of the ^^C properties, and thus the carbon and oxygen production in stellar 
environments, as the function of the N-N interaction strength. The results for 
the resonance energies and widths are shown in Table ^ 

As we mentioned, we correctly reproduced the J — 2]*" energy difference for 
p = 1.0. This required the use of a slightly different force for 2'^ . The resulting 
F^ values should be compared to the experimental figure, 3.7 ±0.5 meVcl. Our 
model performs well, considering the fact that no effective charges were used 
in the calculations. Using an effective charge e^s for both the protons and 
neutrons would lead to F^ multiplied by (1 -I- 2eoff)^. 
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Table 1: The energy, Er (in keV), and gamma width, F-y (in meV), of the -resonance as 
a function of the factor p 



N-N interaction 


MN 


VI 


V2 


P 


e{p) 




e{p) 




e{p) 




1.004 


273.3 


2.85 


294.4 


3.21 


306.0 


3.28 


1.002 


327.5 


2.91 


337.5 


3.26 


343.7 


3.31 


1.001 


353.7 


2.93 


358.7 


3.27 


361.7 


3.32 


1.000 


379.6 


2.96 


379.6 


3.29 


379.6 


3.33 


0.999 


405.2 


2.99 


400.3 


3.31 


397.2 


3.34 


0.998 


430.5 


3.01 


420.8 


3.32 


414.6 


3.34 


0.996 


481.4 


3.07 


460.7 


3.34 


450.0 


3.34 



3 Triple-alpha reaction rate 

The results shown in Table |^ can be used to calculate the triple-alpha reaction 
rates. The reaction rate for the triple-alpha process pmceeding via the ground 
state of ^Be and the Oj-resonance in ^^C is given byEj 

where Ma and Na is the mass and the number density of the a-particle, 
respectively. The temperature of the stellar plasma is given by T. The quantity 
e denotes the difference in energy between the 0^-resonance in -'^^C and the 
3a-particle threshold. The resonance strength ujj is given by 

^7 = -p , p — «r^- (5) 

r Q -r i l ad 

The approximation of the above expression for the decay widths of the 
Oj-resonance follows, because for the a-width Fq,, the radiation width Trad, 
the electromagnetic decay width to the first excited state of ^^C, and the 
electron-positron pair emission decay width Fpair into the ground state of ^^C 
the following approximations hold: (i) Fq ;3> Frad and (ii) Fiad = F^ + Fpair ~ 
F-y. Therefore, Eq. (jj) can be approximated by: 

2Trh^ \ ^ F, 



The two quantities in Eq. (^) that change their value by varying the ef- 
fective N-N interaction are the energy of the 0^-resonance e in -"^^C and its 
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electromagnetic decay width F^. However, the change in the reaction rate by 
varying the effective N-N interaction with e in the exponential factor of Eq. 
is much larger than the linear change with F^. 



4 Stellar burning 

The significance of low-mass (with masses less than about 2 Af©, intermediate- 
mass (between about 2Mq and IOMq), and massive stars (with masses moie 
than about IOMq) in the nucleosynthesis of carbon is still not quite clearllj. 
Some authors claim that low-mass and intermediate-mass stars dominate in 
the production of carbon, whereas others-favor the production of carbon in 
massive stars. In a recent investigation llj by spectral analysis of solar-type 
stars in the Galactic Disk, the results are consistent with carbon production 
in massive stars but inconsistent with a main origin of carbon in low-mass 
stars. The significance of intermediate-mass stars for the production of carbon 
in the Galaxy is still somewhat unclear. Therefore we performed a stellar 
model calculation for a typical low-mass, intermediate-mass, and massive star, 
respectively. 

The calculation of the stellar models are performed with a contemporary 
stellar evolution code, which contains theJatest physics input. In particular, 
using this code, up-to-date solar models Ej can be produced as well as the 
evolution of low-mass stars can be followed through the thermal pulses of 
asymptotic giant branch-(AGB) stars 113. The nuclear network is designed 
preferentially to calculate the hydrogen and helium burning phases in low- 
mass stars. Additionally, the basic reactions of carbon and oxygen burning 
are included, which may destroy the produced carbon and oxygen in massive 
stars. 

In this work the evolution of a 1.3, 5, and 20 solar mass star is calculated, 
which should represent the typical evolution of low-mass, intermediate-mass, 
and massive stars, respectively. The stars are followed from the onset of hy- 
drogen burning until the first thermal pulses in the AGB, or until the core 
temperature reaches 10^ K in the case of the 2OM0 star, as the nuclear net- 
work is not sufficient to go beyond this phase. 

Large portions of the initial mass of a star are returned to the interstellar 
medium (ISM) through stellar winds basically during the thermal pulse phase. 
Unfortunately, basically due to the simple convection model used in stellar 
modeling, the composition of the wind can not yet be determined very accu- 
rately from stellar evolution theory. However, it is beyond the scope of the 
present investigations to determine how and when the material is return back 
to the ISM. Instead we investigate how much C and O is produced altogether 
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Figure 1: The radial abundance profile of carbon and oxygen within a 5.0 shortly before 
the first thermal pulse. The quantities Xc and Xq denote the radial abundances of carbon 
and oxygen, respectively. The models are from left to right with p = 0.996, 0.999, 1.00, 1.001 
and 1.004. 

in the star, which then may be blown away by steUar winds. 

For the three masses quoted above, the evohition of stars is calculated using 
the predictions of the MN interaction with p = 0.996, 0.999, 1.000, 1.001, and 
1.004. In Fig. the radial abundance profiles of C and O are shown for a 5M0 
star, shortly before the first thermal pulse. With the standard cross section 
for the triple-alpha reaction (He-burning front at = 0.95Mq) roughly the 
same amount of carbon and oxygen is produced. Changing the N-N interaction 
strength by only about 0.1%, this ratio could already be altered significantly. 
If the N-N interaction strength is enhanced or reduced by 0.4%, respectively, 
then almost no oxygen or carbon is produced in a 5Mq star. Similar behavior 
can be observed for the 1.3Mq star, while in the case of a 20Mq star, the 
carbon can be destroyed almost totally even for p = 0.999. 

Thus, even with a minimal change of 0.4% in the strength of the N-N 
force, carbon-based life appears to be impossible, since all the stars then would 
produce either almost solely carbon or oxygen, but could not produce both 
elements. For smaller variations in the N-N interaction the development of life 
depends on the stellar population, as for instance a 2OAf0 star with p = 0.999 
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does not produce carbon, while the 1.3Mq and 5Mq stars still do. 
5 Summary 

We have investigated the change of the carbon and oxygen production in the 
helium burning of low-mass, intermediate-mass and massive stars by varying 
the underlying nuclcon-nucleon interaction. A slight variation of the nuclcon- 
nucleon interaction leads to a change of the resonance energy of the second 0+ 
state {O2) in ^^C, thus drastically modifying the of the triple-alpha reaction 
rate. The changes in the relevant nuclear parameters in the triple-alpha reac- 
tion rate, e.g. the resonance energy and radiative decay width of the 0^ state 
in ^^C, were obtained with the help of a microscopic three-cluster model. 

The impact of changing the triple-alpha reaction rate by varying the ef- 
fective N-N interaction in the evolution of 1.3, 5 and 20 solar mass stars is 
calculated by following the typical evolution of a low-mass, intermediate-mass, 
and massive star, respectively. The calculations were carried out with a con- 
temporary stellar evolution code, which contains the latest physics input. The 
result is that even with a small change of 0.4% in the N-N interaction strength, 
carbon-based life appears to be impossible, since all the stars then would pro- 
duce either almost solely carbon or oxygen, but could not produce both ele- 
ments. 

Acknowledgments 

This work was supported by OTKA Grants F019701, F033044, and D32513 
(Hungary), and by the Bolyai Fellowship of the Hungarian Academy of Sci- 
ences. We also acknowledge the support by the Fonds zur wissenschaftlichen 
Forschung in Osterreich, project P10361-PHY. 

References 

1. J. D. Barrow and F. J. Tipler, The Anthropic Cosmological Principle 
(Clarendon Press, Oxford, 1986). 

2. G. K. Opik , Proc. Roy. Irish Acad. A54, 49 (1951). 

3. E. E. Salpeter, Phys. Rev. 88, 547 (1952). 

4. F. Hoyle, D. N. F. Dunbar, W. A. Wenzel, and W. Whaling, Phys. Rev. 
92, 1095 (1953). 

5. C. W. Cook , W. A. Fowler, and T. Lauritsen, Phys. Rev. 107, 508 

(1957). 

6. F. Ajzenberg-Selove, Nucl. Phys. A 490, 1 (1988). 



8 



7. M. Livio, D. HoUowell, A. Weiss, and J. W. Truran, Nature 340, 281 
(1989). 

8. R. Pichler, H. Oberhummer, A. Csoto, and S.A. Moszkowski, Nucl. Phys. 
A 618, 55 (1997). 

9. Y.K. Ho, Phys. Rep. 99, 1 (1983). 

10. A.T. Kruppa, R.G. Lovas, and B. Gyarmati, Phys. Rev. C 37, 383 
(1988); A. Csoto, Phys. Rev. C 49, 3035 (1994); A. Csoto, H. Oberhum- 
mer, and R. Pichler, Phys. Rev. C 53, 1589 (1996). 

11. P.J. Brussaard and P.W.M. Glaudemans, Shell-model applications in nu- 
clear spectroscopy, North-Holland Publishing, Amsterdam 1977. 

12. D.R. Thompson, M. LeMere, and Y.C. Tang, Nucl. Phys. A 268, 53 
(1977); I. Reichstein and Y.C. Tang, Nucl. Phys. A 158, 529 (1970). 

13. A.B. Volkov, Nucl. Phys. 74, 33 (1965). 

14. C. E. Rolfs and W. S. Rodney, Cauldrons in the Cosmos (University of 
Chicago Press, Chicago, 1988). 

15. B. Gustafsson, T. Karlsson, E. Olsson, B. Edvardson, and N. Ryde, As- 
tron. and Astrophys., in press; |astro-ph /98 1 130^ . 

16. H. Schlattl and A. Weiss, Astron. and Astrophys. 347, 272 (1999). 

17. J. Wagenhuber and A. Weiss, Astron. and Astrophys. 286, 121 (1994). 



9 



